Phase Transitions in a Probabilistic Cellular Automaton with Two 
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Abstract 

We study the phase diagram and the critical behavior of a one-dimensional radius-1 two-state 
totalistic probabilistic cellular automaton having two absorbing states. This system exhibits a 
first-order phase transition between the fully occupied state and the empty state, two second- 
order phase transitions between a partially occupied state and either the fully occupied state 
or the empty state, and a second-order damage-spreading phase transition. It is found that all 
the second-order phase transitions have the same critical behavior as the directed percolation 
model. The mean-field approximation gives a rather good qualitative description of all these 
phase transitions. 
PACS: 05.45, 05.90, 03.20 



1 Introduction 

Probabilistic cellular automata (PCA) have been 
widely used to model a variety of systems with local 
interactions in physics, chemistry, biology and social 
sciences [||, |3|, Many of these models ex- 

hibit transcritical bifurcations which may be viewed 
as second-order phase transitions. It has been conjec- 
tured [^, ^ that all second-order phase transitions from 
an "active" state characterized by a nonzero scalar or- 
der parameter to a nondegenerate "absorbing" state 
provided that (i) the interactions are short range and 
translation invariant, (ii) the probability at which the 
transition takes place is strictly positive, and (iii) there 
are no multicritical points belong to the universality 
class of directed percolation (DP). In the Domany- 
Kinzel (DK) cellular automaton ||8|, ^J, which is the 
simplest PCA of this type, Martins et. al. discov- 
ered a new "phase." These authors considered two 
initial random configurations differing only at one site, 
and studied their evolution with identical realizations 
of the stochastic noise. The site at which the initial 



configurations differ may be regarded as a "damage" 
and the question they addressed was whether this dam- 
age will eventually heal or spread. The new phase, they 
called "chaotic" , is characterized by a nonzero density 
of damaged sites. In a recent paper, Grassberger |11] 
(see also Ref. showed that this new phase tran- 
sition also belong to the universality class of directed 
percolation. 



More recently, Hinrichsen |13] stressed that in pres- 
ence of multiple absorbing states with symmetric 
weight, the universality class changes from DP to a 
parity conservation (PC) class, with quite different ex- 
ponents. 

The damage-spreading phase transition can be re- 
lated to usual chaotic properties (i.e. positivity of max- 



imal Lyapunov exponent) of dynamical systems |14|. 
The boundary of the damaged phase depend in gen- 
eral on the algorithm used for the implementation of 
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the evolution rule. Recently, however, Hinrichsen |15] 
formulated an objective, algorithm-independent defi- 
nition of damage spreading transitions. In particular, 
the algorithm with maximal correlation between the 
random numbers used in the updating (i.e. using only 
one random number) allows the computation of the 
minimum boundary of damage spreading transitions: 
inside this region a damage will spread whichever al- 
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Figure 1: Mean- field phase diagram for the density c 
of active sites for p3 = 1, as described by Equation 
(1). The gray levels indicates the asymptotic value of 
c; contours for c = 0, 1/4, 1/2, 3/4 and 1. The dashed 
region marks coexistence of phases. The first order 
phase boundary depends on the initial condition. 



gorithm is used. 

In this paper, we study the phase diagram and the 
critical behavior of a one-dimcnsional radius- 1 totalis- 
tic PCA having two absorbing states, and which ex- 
hibits a multicritical point. In the DK model, which is 
a two-input PCA characterized by the transition prob- 
abilities P(1|00) = 1, P(l|01) = P(l|10) = Pi and 
P(l|ll) = P2) the system has two absorbing states for 
P2 = 1; and exhibit a first-order phase transition at 
pi = 1/2. We may therefore, expect that our system 
will have a nontrivial phase diagram. Moreover, the 
discretization of differential equations lead naturally to 
radius- 1 rules (or three- inputs rules) and we feel that 
it would be useful to have a better knowledge of the 
phenomenology of simple models formulated in terms 
of radius-1 PCA. 



2 The model 



The evolution rule of our model is defined as follows. 
Let s{i,t) denotes the state of the i-th. cell at time t, 
and 9{i, t) = s{i — l,t) + s(i, t) + s(i + 1, t) the sum of 



Figure 2: Mean- field damage-spreading phase diagram 
for p3 = 1. The phase diagram has been obtained 
numerically iterating Equation (2). Levels of damage 
are h = 0, 1/8, 1/4, 3/8 and 1/3. The isolated dot 
marks the multicritical point. 

the cells in the neighborhood, then 



s{i,t + l) 



0, 



if9{i,t) = 0, 
if0{i,t) = 1, 
if9{i,t) = 2, 
ife{i,t) = 3, 



where Xj (j = 1,2,3) is a Bernoulli random variable 
equal to 1 with probability pj, and to with proba- 
bility 1 — Pj- That is, the transition probabilities of 
this model are P(0|000) = 1, P(l|001) = P(l|010) = 
P(l|100) =pi, P(l|011) = P(l|101) = P(l|110) =P2, 
and P(l|lll) = ps- Cell i is said to be "empty" at time 
t if s{i,t) = 0, and "occupied" if s{i,t) = 1. The state 
in which all cells are empty is a fixed point (absorbing 
state) of the dynamics. If J93 = 1. the state in which 
all cells are occupied is also a fixed point. In this work 
we study mainly the case ps = 1, delaying the study 
of the case < 1 to section 5. 

If there exists a stable "active" state such that the 
asymptotic density c of occupied cells is neither zero 
nor one, then the model will exhibit various bifurca- 
tions between these different states as the parameters 
pi and p2 vary. This totalistic PCA can be viewed 
as a simple model of opinion formation. It assumes 
that our own opinion and the opinion of our nearest 
neighbors have equal weights. The role of social pres- 
sure is twofold. If there is homogeneity of opinions, 
then one does not change his mind (absorbing states). 
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Figure 3: Phase diagram for the density of active sites 
c obtained by direct simulations. One run was per- 
formed. Lattice size: 10000, number of time steps: 
10000, resolution: 64 different values for both and 
Vi {P3 = 1)) color codes as in Figure 1. The initial den- 
sity is Co = 0.5. The inset represents the density profile 
along the dashed line. Two critical phase transitions 
are evident. 



otherwise one can agree or disagree with the majority 
with a certain probability. 

Let us assume that one starts with a random con- 
figuration half filled with O's and I's. On the line 
Pi + P2 = 1) if Pi is small, to first order in pi, clusters 
of ones perform symmetric random walks, and have 
equal probabilities to shrink or to grow. Slightly off 
this line, the random walks are no more symmetric, 
and according to whether pi is greater or less than 
1 — P2, all the cells will eventually be either occupied 
or empty. pi+P2 = 1 is, therefore, a first-order transi- 
tion line. This line, however, cannot extend to pi = 1 
since, for p2 = 0, our model is similar to the diluted 
XOR rule (Rule 90) of the DK model. We thus expect 
that our model will exhibit, for a certain critical value 
of pi, a second-order phase transition coinciding with 
the damage-spreading phase transition. For pi = 1 
and p2 = 0, frustration is large, and our model is just 
the modulo-2 rule (Rule 150) as in the DK model. 

For p3 = 1, the model is symmetric under the ex- 
change pi 1 — p2, and <-> 1. Thus, along the line 
P2 = 1 — pi the two absorbing states and 1 have 
equal weight, and, according to Hinrichsen |l^, this 
transition should belong to the PC universality class. 



Figure 4: Profile of the hysteresis region for several 
values of the noise e and relaxation time T. Data ob- 
tained from the local structure approximation of order 
I = 6; the regions corresponds to the intersections of 
the hysteresis cycle with c = 0.5 (horizontal dashed line 
in the inset). The hysteresis boundary lines smoothly 
join at pi = 0, p2 = 1 (not represented). Darker 
to lighter areas correspond to T = 500, e = 0.0001; 
T = 1000, £ = 0.0001; T = 500, e = 0.001. The 
dashed lines represent the mean-field hysteresis region. 
The inset represents the cycle along a line parallel to 
the diagonal pi = P2- 

3 Mean-field approximation 

In order to have a qualitative idea of its behavior, we 
first study our model within the mean-field approxi- 
mation. If c{t) denotes the density of occupied cells at 
time t, we have 



c{t + I) = ^pic{t) {I - c{t)Y + 

2,p2C^{t){l-c{t)) + c\t). 

This one-dimensional map has three fixed points: 

3pi-l 



(1) 



0, 1, and c* 



1 + 3pi - 2>p2 



is stable if pi < ^, 1 is stable if p2 > |j and c* is 
stable if > ^ and p2 < §• In the (pi,p2)-parameter 
space, the bifurcations along the lines 



1 2 

Pi = 3, 0<p2<2, 



1 2 

3<Pi<l, P2 = 3, 
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Figure 5: Phase diagram for the damage spreading 
from direct numerical simulations. Numerical values 
as in Figure 3, color codes as in Figure 2. Traces of 
the second order phase transitions are present; they 
join to the first order (cq = 0.5) phase boundary. 

may be viewed as second-order transition lines; the 
first one between the active and empty states, and 
the second one between the active and fully occupied 
states. 

In the domain defined by 

1 2 

< pi < - and - < P2 < 1 

fixed points and 1 are both stable. Their basins 
of attraction are, respectively, the semi-open intervals 
[0, c*[ and ]c*,l]. Therefore, if we start from a uni- 
formly distributed random value of c, as time t goes to 
infinity, c{t) tends to with probability c*, and to 1 
with probability 1 — c*. Since, for pi +p2 = l, c* = |, 
the line defined by 

1 2 

Pi+P2 = 1, with < pi < - and - < P2 < 1 

o o 

is similar to a first-order transition line, although we 
cannot define a free energy for our dynamical system. 

In general, first order phase transitions are associ- 
ated to the presence of an hysteresis cycle, due to the 
coexistence of two phases, both in equilibrium. The 
problem of a proper definition of such a cycle is rather 
subtle: since this model presents absorbing states, it 
is out of equilibrium and the ergodicity is always bro- 
ken, even for finite-size systems. Thus, once the sys- 
tem settles into one of the two absorbing states, it 
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Figure 6: Mean-field phase diagram for = 0. Color 
codes as in Figure 1. 

never exits, even if it is not stable. One can circum- 
vent this problem adding a little of noise, i.e. setting 
P(0|000) = 1 - P(l|lll) = £. This assumption, how- 
ever, brings the model into the class of equilibrium 
models, thus forbidding the presence of a true phase 
transition in one dimension. However, if e is small, 
the system can get trapped for a long time T in a 
metastable region near an absorbing state. The time 
T is also a function of the system size L. Thus, we 
have to perform carefully the limits L ^ oo, T ^ oo 
and e — > 0. In practice, we have observed that the 
intensity of noise e strongly affects the amplitude of 
the hysteresis region, while, for a system size L suffi- 
ciently large, there exists a large interval of possible 
time length T. The results from numerical simulations 
are reported in section 4. 

In the mean-field approximation the hysteresis re- 
gion ranges from < < 1/3, p2 = 2/3 to pi = 1/3, 
2/3 < P2 < 1- The mean-field phase diagram for the 
density is shown in Figure 1. 

One can write down the mean-field equation for the 
damage spreading taking into consideration all pos- 
sible local configurations of the two lattices. Let us 
denote by ri{t) the density of damaged sites at time t. 
The evolution equation for rj depends on correlations 
among sites, i.e. on the order of the mean field for the 
density. Using the simplest factorization for the den- 
sity, r]{t) depends on c{t). The evolution equation for 
the minimum damage ij, i.e. the damage when the evo- 
lution of the two replicas is computed using only one 
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Figure 7: Mean-field damage-spreading phase dia- 
gram, for p3 = 0. The diagram has been obtained 
numerically iterating Equation (2). 

random number, is given by 
ri{t + l)= ^ 7r(c(t),siS2S3)7r(r?(t),/ii/i2/i3) • 

hih2h3 l^J 
|-P(1|S1S2S3) - -P(0|S1S2S3 © hih2h^)\ 

where 

7r(a,xiX2X3) = a^i+^2+^3(l - 

and the symbol © represents the bitwise sum modulus 
two (XOR) of two Boolean configurations. The value 
for c{t) is given by mapping (|^). 

We have numerically iterated Equation (^). The 
resulting mean-field phase diagram for the damage is 
shown in Figure 2. 

4 Numerical simulations 

We used the fragment method [^] to determine the 
phase diagram. Figure 3 represent a plot for the 
asymptotic density c of occupied cells when the ini- 
tial fraction of the occupied cells is equal to 0.5. 
The scenario is qualitatively the same as predicted by 
the mean-field analysis. In the vicinity of the point 
{PI1P2) = (0, 1) we observe a discontinuous transition 
from c = to c = 1, while close to {pi,P2) = (1, 0) the 
transition is continuous. The two second-order phase- 
transition lines from the non fully occupied states to 
either the fully occupied state or the empty state are 
symmetric and the critical behavior of the respective 
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Figure 8: Cut of the density phase space along the pi- 
P2 plane for different values of p3. From left to right 
and for bottom to top: p3 = 0.0, 0.2, 0.4, 0.6, 0.8 and 
1.0. Color codes as in Figure 1 

order parameters, 1 — c and c, are the same. We have 
checked that the critical exponent f3 for these transi- 
tions far from the crossing point are numerically the 
same of DP. 

These two transition lines meet on the diagonal 
Pi + P2 = 1 and become a first-order transition line. 
Crossing the phase boundaries on a line parallel to the 
diagonal pi = P2, the density c exhibits two critical 
transitions, as shown in the inset of Figure 3. Ap- 
proaching the crossing point, the critical regions of the 
two transitions vanishes, the transition itself become 
sharper and the corrections to scaling increase. 

At the crossing point, around (0.5,0.5), the two at- 
tractors have symmetrical weight. If we relabel couples 
of I's with the symbol 1, couples of O's with the symbol 
and couples 01 or 01 with the symbol A, we fulfill the 



condition stated by Hinrichsen |13| to have symmetric 
absorbing states, whose cluster are always separated 
by an active {A) layer. We have performed prelimi- 
nary measurements of the exponent (3 for the density 
of A couples along the line pi + P2 = 1, and obtained 
(3 ~ 0.65(5), which is consistent with the values for the 
PC process reported in Refs. |l^, |l^]. The transition 
point is at pi = 0.460(2), which defines the position of 
the tricritical point. 

In order to study the presence of an hysteresis re- 
gion, we performed several scanning of the first-order 
transition using the local structure approximation |19|. 
We cut the phase diagram with a line parallel to the 
diagonal pi = P2, increasing the value of pi and p2 af- 
ter a given relaxation time up to p2 = 1; then reverting 
the scanning up to pi = 0. The results are reported in 
Figure 4. The size of the hysteresis region grows with 
the level of the noise e. Preliminary simulations show 
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that the relaxation time T and the noise level e scales 
as eT^'^^^ = const. 

The phase diagram for the damage spreading is 
shown in Figure 5. Outside the true damage-spreading 
region, there appear small damaged domains on the 
other phase boundaries. This is due either to the di- 
vergence of the relaxation time ( second-order tran- 
sitions) or to the fact that a small difference in the 
initial configuration can drive the system to a different 
absorbing state (first-order transitions). 

The "chaotic" domain near the point {pi,P2) = (IjO) 
is stable regardless of the initial density. Note that on 
the line P2 = its boundary coincides with the transi- 
tion line from the active state to the empty state. Us- 
ing an argument similar to the one in Ref. [1T2|| would 
prove that the derivatives of the two boundary curves 
coincide, that is, the chaotic phase exhibits a reentrant 
behavior. Here again, the symmetry of the phase dia- 
gram implies a similar behavior on the line pi = 1. 

5 Extended phase diagram 

It is interesting to study our model when the transition 
probability P(l|lll) = p^ < 1. In this case, there is 
only one absorbing state. 

5.1 Mean-field approximation 

For some values of p^, the mean-field phase diagram 
still exhibits a first-order phase transition. Let us first 
examine in detail the particular case p^ = 0. As for 
P3 = 1, the solution c = looses its stability at p = 
1/3. The other solutions for the density c are the roots 



Figure 10: Cut of the density phase space along the 
Pi-ps plane for p2 = 1. Color codes as in Figure 1. 

of the quadratic equation 

c^{pi -p2) + c{p2 - 2pi) +pi - i = 0. 

These two roots are real if + 4{pi — P2)/3 > 0. This 
is always the case if pi > 1/3, while for pi < 1/3, p2 
must lie outside the interval [|(1 — ^/T^^3pl), |(1 + 
^/T^^3pl)~\ . Since the density c has to be such that 
< c < 1, only the condition I > P2 > |(l + \/l ~ 3pi) 
is meaningful. These solutions are stable. Therefore, 
there is a domain in which two stable solutions co- 
exist. The boundaries of this domain are the three 
straight lines pi = ^, p2 = I andp2 = |(l + -v/l — 3pi), 
where for this last line p2 belongs to the interval 
[|,^]. When, inside the domain, we approach the 
point (pi,P2) = (^jf); the density c tends contin- 
uously to zero, which shows that, at this point, the 
transition is second-order. The corresponding phase 
diagram is represented in Figure 6. 

For the chaotic phase, the boundary of the domain 
of stability is shown in Figure 7. Close to the point 
(Pi)P2) = (liO) the boundary is almost the same as 
for p3 = 1 but there exists another domain close to the 
point {pi,P2) = (1, 1) which did not exist for p^ = 1. 
In the vicinity of (1, 1), very small clusters, which tend 
to grow, are then fragmented because ps = 0. We have 
also studied the case < ps < 1. 



5.2 Numerical simulations 

Our numerical simulations give results in good qual- 
itative agreement with the mean-field approximation. 
Figure 8 shows different phase diagrams in the {pi,P2)- 
plane corresponding to different values of p^. As ex- 
pected, the boundary of the active phase (c 7^ 0) for 
P2 = does not vary much with p^ since, in this case, 
the probability to find a cluster of three or more one 
is vanishing. No first order phase transition for the 
density is apparent by these simulations. 

Figure 9 shows the domain of stability of the chaotic 
phase in the (pi,p2) -plane for different values of p^. 
The domain of stability of the chaotic phase in the 
(pi,P3)-plane for p2 = 1 is represented in Figure 10. 
One can see that there are damage-spreading domains 
corresponding to the density phase boundaries. 

6 Conclusion 

We have studied a radius-1 totalistic PCA whose 
transition probabilities are given by P(1|000) = 0, 
P(l|001) = P(l|010) = P(l|100) = pi, P(l|011) = 
P(l|101) = P(l|110) =p2, and P(l|lll) = Ps, within 
the framework of the mean-field approximation and 
using numerical simulations. 

For p3 = 1, the system has two absorbing states and 
four different phases: an empty state in which no cell 
is occupied, a fully occupied state, an active phase in 
which only a fraction of the cells are occupied, and a 
chaotic phase. This system exhibits: 

1. A first-order phase transition from the fully occu- 
pied state to the empty state. 

2. Two second-order phase transitions from the ac- 
tive state to the fully occupied state and from the 
active state to the empty state. 

3. A second-order phase transition from the active 
state to the chaotic state. 

For p3 < 1, our numerical simulations show that 
there is no first-order phase transition while the mean- 
field approximation still predicts the existence of such 
a transition. 
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